


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1959 


A graphical-numerical prediction of the 
1000-mb surface with non-adiabatic warming. 


Wang, Yeh-chun. 


Monterey, California: U.S. Naval Postgraduate School 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


1 ih D U DLEY research materials and institutional publications created by the NPS community. 
۱ Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
d 
KNOX appointed — and published — scholarly author. 


WU LIBRARY Dudley Knox Library / Naval Postgraduate School 


411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 








http://www.nps.edu/library 


NPS ARCHIVE 
1959 


WANG, Y. 





A GRAÞHICAL-NUNERICAL PREDICTION. 
OF THE 1000-MB SURFACE WITH M 


NON - ADIABATIC WARMING ge 


тил "ray 


. fEH-CHUN WARG 


Bê 


LIPRARY 
U.S. NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 

















A GRAPHICAL-NUMERICAL PREDICTION OF THE 


1000-MB SURFACE WITH NON-ADIABATIC WARMING 


e? 


Yeh-Chun Wang 








A GRAPHICAL-NUMERICAL PREDICTION OF THE 


1000-MB SURFACE WITH NON-ADIABATIC WARMING 


by 


Yeh-Chun Wang 
” 


Lieutenant Commander 


The Republic of China Navy 


Submitted in partial fulfillment of 
the requirements for the degree of 


MASTER OF SCIENCE 
IN 


METEOROLOGY 


United States Naval Postgraduate School 
Monterey, California 
United States of America 


DYS 








DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTBREY CA 93943-5101 


A GRAPHIC‘ L-NUMERICAL PREDICTION OF THE 
1000-MB SURFACE WITH NON-ADIABATIC WARMING 
by 


Yeh-Chun Wang 


This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 
METEOROLOGY 
from the 


United States Naval Postgraduate School 














ABSTRACT 
A two-level graphical-numerical prediction model is extended so as 
to incorporate some orographical and non-adiabatic warming or an 
effects in a manner such that almost no extra work is added. 
The model is applied to several cases where non-adiabatic heat 


exchanges are pronounced. The resulting prognoses show a material improve 


ment over a similar model which omits the non-adiabatic influence. 
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1. Introduction 

It is a -ell-kno-n fact that during the winter season cold cP air 
flowing out from continent over the warm ocean surface will acquire heat 
energy from below so that the lower layers of the atmosphere are warmed. 
As Craddock [1] measured in the case when the arctic air flowed out from 
Iceland to the British Isles, the heating rate may be as large as 65 /1y 
per hour and in the course of a day the 1000-500 mb layer thickness could 
increase approximately 800 feet. Hence, in such cages, an appreciable 
error will be produced in the results of those numerical prediction models 
which assume that the air motion is adiabatic. Several techniques have 
been introduced to eliminate this non-adiabatic error; Reed [5] has 
incorporated Burke's semi-empirical treatment for the air mass transforma- 
tion cP-mP into a graphical model of numerical prognosis. His results show 
considerable improvement over a model omitting the non-adiabatic heating. 

In the present investigation, the non-adiabatic heating (or cooling) 
is approximated in a rather crude but extremely simple manner. As a 
result there is no additional work when non-adiabatic heating and the 


orographic influences [3] are included. 





Development of the prediction equation 


The basic vorticity equation in the (x,y,p,t) coordinates takes the 





form 


DF +t $2 _ E (1) 
им A 4 


where je is the individual rate-of-change operator following the isobaric 
projection of the three-dimensional motion, 
J the relatíve vorticity, 
£ coriolis force, 
WV del operator. 
Both | and Y refer to values measured on an isobaric surface. 
In Eq. (1) the terms involving the vortex tube and vertical advection 
of vorticity have been neglected, and also the term 7: bel: 


Expanding the left side of Eq. (1), 
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By equation of continuity, 


= 6 
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Eq. (2) becomes 


Duvet) 5 3q (3)‏ اج 


where 
Ap f 
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zê WAR 
Assume now that w has a vertical parabolic distribution plus the 


terrain-induced vertical motion. 
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Here : W are the surface pressure, density and wind 
Kei M , 
velocity; H is the height of the terrain; /ris an upper-pressure level 
$ 
(to be taken as 500 mb) where vertical motion is strongest; and 2 is a 
constant of proportionality. 
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The term on the right side of Eq. (4) involving Y -YH represents the 
terrain induced vertical motion. Here the surface vertical velocity A, is 


assumed to be proportional to the horizontal wind and the gradient of 


terrain height H at the surface; i.e. ZŁA y H . Moreover, 


since girade E: it follows that 


A), = — es لا‎ YH 


Finally, the terrain-generated vertical velocity is assumed to 
decrease linearly with pressure, and thus we arrive at the first term consti- 
tuting Jin Eq. (4). 

The second term of Eq. (4) represents the commonly used parabolic dis- 
tribution of vertical velocity assumed to be typical of large-scale press- 
ure systems [4]. Where terrain effects are negligible, this term corre- 


tal 


sponds to low-level ''isobaric ' velocitv divergence and high level con- 
vergence (or vice versa), together with the level of non-divergence at A 


Differentiating Eq. (4) with respect to ZZ, substituting it into 


Eq. (3) and evaluating the result at "ды gives 
о 
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The next step consists of utilizing the first law of thermodynamics 
to provide a second equation involving Cn: The non-adiabatic heat 
addition to the surface-500 mb laver will be assumed to be proportional 
to the rate at which the surface air parcels are exposed to a changing 
temperature of the underlving surface. 

Expanding the individual change of the potential temperature, -) ; 


and dividing through by 2 : 
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The subscript p at the W operator will be omitted henceforth. 


It can be shown from Poisson's equation 
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that with pressure held constant 
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where e is the specific volume and R, the universal gas constant. Sub- 


stituting these identities into Eq. (6) yields 
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From the hydrostatic equation, 
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substitute into Eq. (7) after multiplying through by A» giving 
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It follows that 
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Equation (8) states that the local change of the rate of change of the 
pressure height with respect to presgure is the result of: a) advection, 
b) vertical motion, and c) the non-adiabatic heating process represented 
by the last term. 

At this stage, the non-adiabatic addition of heat will be assumed 
proportional to the expression V 7 where VA is the wind at the 
surface and T. is the surface temperature of underlying ocean surface. 


Hence Eq. (8) may then be expressed in the form 


9 /, 
GE EE = Kw. vr (9) 


— is the constant of proportionality. Strictly speaking, f. is not 
a constant but may possibly vary with synoptic situations and also vary 
Slowly with time. However, £ vyill be assumed constant for the present 
investigation. 


Nov substitute Eq, (4) into Eq. (9) and make geostrophic approximation in 


the advection term. 
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By integrating between ji and ft assuming O to be a constant, we obtain 
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P (Aw. VPI) dp. 


Here h is the thickness between and 22 
[rs pi 


In the last term of the preceding equation, W and T, are independent of 


pressure, hence 


| 
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where / is a new constant. 


Solving Eq. (10) for GL, 


and substituting it into Eq. (5) gives 
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Rearranging, 
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Nov, we make the geostrophic approximation for Bi in Has (il). 
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in comparison to (2-2). With multiplication by APA | e obtaín 
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The coriolis force term on the right side of Eq. (12) has been 


manipulated as follows: 
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Therefore Eq. (12) becomes 





= O 3 -2,+Ch+G +F+NT] (13) 


The function// [varies very slowly with time and may be assumed 
to be a constant for the normal forecast period. The functions G and 
F depend on geographical position only and are invariant with time, 
hence they may be inserted on the left side of Eq. (13). Equation (13) 
shows that the quantity ( Z-Z, + J + G + E + VT.) is con- 
servative with respect to wind field Se By using the geostrophic 


о 
approximation for\Y , i.e. 
Ó 


E 2 l _ 7 
grad ij 6778 


Eq. (13) mav be written in Jacobian form as 
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In order to arrive at a more conservative field for advecting the 
conservative property, the Jacobian identity, J(x,y) = J(x+y,y), may be 


applied to Eq. (14), giving 


+F ANT]‏ +۰ 2-2 اد 
FI 12۰64 Eent, Set) 09‏ == 


The use of this form implies that the / Z, - CZ FG FF 4 WT; ] field 
is expected to be more conservative than the CS field for advection pur- 
poses. This is not necessarily always the case since the thickness field 
appears in the advecting quantity. Nevertheless, in the examples tested, 


Eq. (15) gave good results. 





To depict the local change of [ Z,-2. TONS T F-t^T ] ve 
shall advect this function for a short time interval, say twelve hours 
or twenty-four hours. Denote this change by A . Since G, F and NT are 


not functions of time, 
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Taking the space average over Eq. (17) and substituting back to Eq. (17), 


vields 
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Averaging n times, ne 
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47 др ——— 
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However, this series converges rapidly; and, as shown by Fjortoft [2], 
it can be approximated as 
, Dm 
MA TE S | 
a TOBI 12 (19) 
In the cases tested in this investigation, it was found that the 
second term on the right side of Eq. (19) can be neglected except when 
a marked change is indicated. In the routine work where a time saving is 
of great importance, AZ = 6 can be used without much sacrifice in 
accuracy. In utilizing the high speed electronic computer, Eq. (16) 


a 


could be solved to anv desired accuracv. 
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3. Procedure 

In order to provide a comparison to Reed's results [5], the same 
case has been chosen for prognostication. Two other series were also 
chosen for testing purposes. One series has the storm located in the 
Gulf of Alaska. Each series contains three consecutive periods with a 
time interval of twelve hours. Except in the same case as Reed's, where 
the prognosis without heating effect was already available, two prog- 
noses for each case have been obtained, one omitting and the other includ- 
ing the heating effect. 


By omitting the heating effect, Eq. (15) becomes 
2 27 un ۱ m. 
die Zt Cf Eres FE ] 


=- FILE +A GA E-E Ce Ge ] cn 


The following steps are carried out in each prognostication: 

(1). A space-mean chart for 20 is graphically constructed using a 
grid distance of 600 km. 

(2). With C = 0.5, the current thickness chart is multiplied by 
1/2 (this is achieved by simply tracing every other line on the thickness 
chart) and added to z. 

(3). Next the result of step (2) is added to the (G + F) function 


provided by Haltiner and Hesse [3]. 
loo JT 

© 

A 


land area of Alaska (see Section 4), a (NT. ) function chart was constructed 


over the ocean and N = -0.05 over the 





(4). With N ==0 / 


based on the monthly mean ocean surface isotherms published by the Hydro- 
graphic Office and the current reported temperatures over the land. 
(Fig. 19). This chart of ؟‎ ) 1 „as added to the result of (3). Note 


this step is unnecessary if non-adiabatic heating is omitted. Also it is 
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suggested, since monthly mean surface temperature isotherms over the 
ocean are used, a chart oí (G + F 4 NT. ocean) function may be prepared 
for each winter month. This will combine step (4) into step (3) if the 
(NT) function over land is also approximated by the monthly mean iso- 
therms over the land. Up to this step, a chart or function of 
(2, + Ch + G + F + NT |) has been obtained and will be referred to as 
steering chart henceforth. 
(5). The 2o chart is subtracted from steering chart to obtain the 
quantity to be advected. 
(6). Next advect the quantity obtained in step (5) with the steering 
chart for twelve-hour intervals. 
(7). Subtract the result of step (6) from the result of step (5) to 
get A, 
(8). Add 1/242, to A and multiply the result by 2/3 to obtain B. 
(9). Average B to get B and multiply by 4/3. 
(10). Add the result of (9) to (8), to get AZ. 
(11). AddA 20 to the initial 1000-mb chart to give the prognostic 


1000-mb chart of twelve hours later. 
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4. Discussion 
In this investigation, three constants are of major importance, 
namely; C, k, and N. Otr these three constants, k, has been discussed in 
the Haltiner and Hesse paper [3], and only C and N will be discussed here. 
(1). The constant C. 
several values of C have been used by various investigators. 
In the present investigation, the value of C = 0.5 fits best in the cases 
tested. Values lower than this will give an advection field too weak, 
While higher values give an advecting field too strong. 
(2). The constant N. > 
In this investigation the primary interest is in the area of 
outflow of cP air from inland Alaska to the sea. Over the ocean the value 


13041 


= In the non- 


of constant N was found empirically to be -0.1 
adiabatic modification of an air mass, the processes which actually take 
place are not simple. The degree and rate of the modification are depend- 
ent upon many factors such as: a) the temperature difference between the 
cP air mass and sea surface; b) the roughness of the sea surface; c) the 
wind speed over the sea surface; d) the length of time the air mass has 
been over the ocean; e) the stability of the air mass; f) the moisture com 
tent of the air mass; ete. The inter-relations among these factors are not 
well known quantitatively; consequently, it is difficult to determine in a 
specific case the relative importance of the various factors. For example, 
air moving with a higher speed is not necessarily being modified pro- 
portionately more than air moving at a much lower speed. The reason for 
this is that the modification is not only a function of the distance over 
which the air has traveled from the coast, but also a function of the time 
the air mass has been over the ocean. The expression Zeen implies 
that when the value of Y is doubled, the amount of modification will 
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double, which is not necessarily true. Thus, the determination of the 
value of N is a matter of compromise of the many different influences. 
For the relatively small sample tested, the value of N = -0.1 appears to 
be the best choice thus far. However, application of this value shows 
considerable improvement over the other models which treat the non- 
adiabatic heating effect in a different way. 

In the course of this investigation, it appeared that the advecting 
field tends to move the high pressure cell down stream over inland 
Alaska, with the result that the prognostic pressure height was too high 
over the southwestern side of the High and too low over the southeastern 
side of the High. For example, in the case 1500 GMT January 4, 1956 to 
0300 GMT January 5, the high was built up on the east side. In this situa- 
tion, an attempt was made to apply the non-adiabatic effect over inland 
Alaska by putting N = -0.05, obtaining a £(NT,) based on the isotherms on 
the current surface map. Considerable improvement resulted. Such an 
approach provides a continuous temperature field from the inland out to the 
sea. This value of N was also applicable to the other inland cases. The 
combined functions (NT) for the land and the sea are shown in Figs. 19, 
20, 21, 22. Frem these figures, it seems feasible to use the monthly mean 
temperature over Alaska land area as the basis of f(NT,) over land, so that 
step (4) can be entirely incorporated into step (3) by preparing a single 
chart ot f(G + F + NT.) for each winter month. Thus a considerable amount 
of time can be saved in daily routine work. Since the non-adiabatic cool- 
ing process over the land, probably mostly due to radiational cooling, is 
beyond the scope of this investigation, no further discussion will be given. 

In addition to the empirical determination of constants, another 
aspect at the advecting stage in the procedure warrants some discussion. 


As the air mass moves southward over an increasing temperature field, the 
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4, 





air mass will be warmed up from below. An assumption was made that the 
air parcel directly in contact with the sea surface will immediately 
acquire the same temperature as that of the sea surface. In this case, 
the air mass tends to be unstable. This instability will enhance the 
upward heat transport, and so increase the modification. On the contrary, 
when the air mass moves northward over a decreasing temperature field, 

by the same assumption it tends to be stabilized. In the first case, the 
thickness of the lower layer will increase and in the latter case the 
thickness will decrease. However the magnitude in the first case of 
increasing thickness is normally much larger than the magnitude in the 
latter case of decreasing thickness. Because of this consideration, an 
attempt was made to advect less distance in cold advection and to advect 
the full distance in warm advection. However, this did not improve the 
results, since the property to be advected is a combination of 20. Za 

Ch, G, F, and NT > and not a function of NT alone. As a matter of fact, 


since the value of N has been determined empirically, it represents a 


compromise of all the various influences. 
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Results 





The results of the tests are displayed in Table I which shows the com- 
parison of the series, with and without non-adiabatic warming or cooling. 
The statistical parameters computed were based on a seventy-point grid of 
which forty-seven points were over the sea. These grid points are the 
intersections of the meridians and latitude parallels which are 5° apart 
extending from longitude 165° 4 to longitude 120°W and from latitude 65°N 
down to latitude 30°N. Due to the convergence of meridians, the ed 
are denser over the northern part of the area, Since the area of denser 
points is inland and coastal Alaska, which is a more critical area, the 
values in Table 1 are more conservative and do not exaggerate the improve- 


ment due to the inclusion of the non-adiabatic heating. 
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Figs. 1 to 16 indicate the results of four casts. These are: 
Case 1) 1500 GMT January 4, 1956 to 0300 GMT January 5, 1956; Cases 2, 
3, 4) 0000 GMT January 13, 1959 to 1200 GMT January 14, 1959, consisting 
of three consecutive prognoses. The other three cases from 0000 GMT 
January 13, 1959 to 1200 GMT January 14, 1959 are not shown. It has also 
been tiwi that the addition of term NT, gives no improvement when a 
cyclone comes from the central North Pacific area; however, neither is 
there a decrease in accuracy when f(NT,) is included, 

It can be seen from the figures that the greatest improvement in the 
prognostic pressure heights is found over the area on the west of the storm 
which is the normal location of the cP air. This was true in all the cases 
tested. 

On the basis of the sample tested it may be concluded that this model 
shows very good promise. However, further testing would be desirable to 
determine the best values of the various constants. The non-adiabatic 
heating term can readily be incorporated into the programming of baroclinic 
models for the electronic computer on which further testing would be 
greatly facilitated. There is otwiously a considerable saving in time 
compared to the Reed method when this model is applied manually in the 


graphical technique. 
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